Tune-out wavelengths for potassium 
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The five longest tune-out wavelengths for the potassium atom are determined using a relativistic 
structure model which treats the atom as consisting of a single valence electron moving outside a 
closed shell core. The importance of various terms in the dynamic polarizability in the vicinity of 
the 4pj, 5pj and 6pj transitions are discussed. 
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I. INTRODUCTION 

The dynamic polarizability of an atom gives a measure 
of the energy shift of the atom when it is exposed to an 
electromagnetic field [l|, For an atom in any given 
state, one can write 
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where ^^(w) is the polarizability of the quantum state 
at frequency uj, and is a measure of the strength of 
the AC electromagnetic field. The limiting value of the 
dynamic polarizability in the a; — >■ limit is the static 
dipole polarizability. 

The dynamic polarizability will go to zero for certain 
frequencies of the applied electromagnetic field. The 
wavelengths at which the polarizability goes to zero are 
called the tune-out wavelengths [J Atoms trapped 
in an optical lattice can be released by changing the 
wavelength of the trapping laser to that of the tune- 
out wavelength for that atom. Very recently, tune-out 
wavelengths have been measured for the rubidium and 
the potassium atoms d, 0. The advantage of a tune- 
out wavelength measurement is that it effectively a null 
experiment, it measures the frequency at which the po- 
larizability is equal to zero. Therefore it does not rely 
on a precise determination of the strength of an electric 
field or the intensity of a laser field. 

In the present manuscript a calculation of the five 
longest tune-out wavelengths for the potassium atom is 
presented. The method applied is a fully relativistic ver- 
sion of a non-relativistic semi-empirical fixed core ap- 
proach that has been successfully applied to the descrip- 
tion of many one and two electron atoms 

An 

extensive discussion is made about those parts of the os- 
cillator strength sum rule that have the largest influence 
in the determination of the tune-out wavelengths. 



II. FORMULATION 

The calculation methodology is as follows. The first 
step involves a Dirac-Fock (DF) calculation of the K+ 



TABLE I: The cutoflF parameters, 
potential. 



of the core polarization 



e 


J 


pej (ao) 


s 


1/2 


2.1360 


p 


1/2 


2.0324 




3/2 


2.0289 


d 


3/2 


2.3610 




5/2 


2.3633 



ground state. The single particle orbitals are written as 
linear combinations of analytic basis functions. The func- 
tions chosen are the S'-spinors introduced by Grant and 
Quiney [III 113 ■ S-spinors can be regarded as relativistic 
generalizations of the familiar Slater type orbital. 

The effective interaction of the valence electrons with 
the core is then written 



H = ca ■ p + iSmc^ + Vcorc{r) 



(2) 



where m is the electron mass, c is the speed of light, p is 
the momentum operator, a and (3 are 4x4 matrices of 
the Dirac operators The core operator is 



K:oro(r) 



Z 



+ ydh-(r) + Kxc(r) + yp(r) . (3) 



The direct and exchange interactions of the valence elec- 
tron with the DF core were calculated exactly. The £- 
dependent polarization potential, Vp, was semi-empirical 
in nature with the functional form 
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The factor, acore is the static dipole polarizability of 
the core and gjj{r) = 1 — exp(— r^/p^^) is a cutoff func- 
tion designed to make the polarization potential finite at 
the origin. The cutoff parameters, pgj were tuned to re- 
produce the binding energies of the ns ground state and 
the np, and nd excited states and are listed in Table HI 

The effective Hamiltonian for the valence electron was 
diagonalized in a large L-spinor basis L-spinors 
can be regarded as a relativistic generalization of the 



2 



Laguerre type orbitals that are often used when solv- 
ing the Schrodinger equation This basis can be en- 
larged towards completeness without any linear depen- 
dence problems occurring. There is effectively no er- 
ror due to the incompleteness of the basis set in the 
present calculation. The present relativistic configura- 
tion interaction plus core polarization calculations typi- 
cally used 50 positive energy and 50 negative energy L- 
spinors for each (£, j) symmetry. This approach is named 
the relativistic configuration interaction plus core polar- 
ization approach (RCICP). The non- relativistic approach 
from which the method is derived is called the configu- 
ration interaction plus core polarization (CICP) method 
For the purpose of comparison, we present results 
of calculations using the all-order single-double imple- 
mentation of relativistic many body perturbation theory 
[IKil] (MBPT-SD). The area of commonality between 
the RCICP and MBPT-SD approaches is that both cal- 
culations have minimal numerical uncertainties. The two 
methods use different approximations to treat the inter- 
action with the core and core-valence correlations, but 
the subsequent calculations within their respective the- 
oretical frameworks have effectively no significant errors 
due to basis set incompleteness. 

III. RESULTS 

A. Energies 

Table [IT] gives the energies of some of the low lying 
states of potassium. One of the interesting aspects of the 
table concerns the spin-orbit splitting of the 5pj and 6pj 
states. The polarization potential parameters Pi.1/2 and 
Pi, 3/2 were tuned to give the correct spin-orbit splitting 
of the 4pj states. Making this choice resulted in the 
spin-orbit splittings for the bpj and (Spj states also being 
very close to experiment. Similarly, tuning the P2,d./2 
and P3^5/2 parameters to give the correct 3dj spin-orbit 
splitting also resulted in the spin-orbit splittings for the 
Adj and bdj levels also in agreement with experiment. 

B. Dipole matrix elements 

Table imi gives the reduced matrix elements for a num- 
ber of the low-lying transitions of the potassium atom. 
These matrix elements were computed with a modified 
transition operator [1,[H,[i3j e.g. 

r = r-(l-exp(-rVp^))'^'^ (5) 

The cutoff parameter used in Eq. ([S]) was 2.176 oq, the 
average of the s, p and d cutoff parameters (note, the 
weighting of the s was doubled to give it same weighting 
as the two p and d orbitals). These matrix elements 
are compared to the MBPT-SD matrix elements [l3llT8|. 
Some reduced matrix elements derived from experiment 



TABLE II: Theoretical and experimental energy levels (in 
Hartree) for some of the low-lying states of K. The energies 
are given relative to the energy of the K"*" core. The experi- 
mental data were taken from the National Institute of Science 
and Technology (NIST) tabulation 
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are also given |13|. It should be noted that there are 
some small differences between the most recent MBPT- 
SD reduced matrix elements [l3l | and earlier work using 
this method flQ^ . 

The agreement between the RCICP and MBPT-SD 
calculations for the larger matrix elements is very good. 
The two calculations agree to better that 1% for the 
4s — 4p, 4p — 5s, 5p — 6s and Sd — Ap transition man- 
ifolds. The same level of agreement is achieved for the 
non-relativistic CICP matrix elements. The CICP ma- 
trix elements are taken from earlier calculations of dis- 
persion coefhcients for atomic pair involving potassium 
11 IMm. The agreement between the CICP, RCICP 
and MBPT-SD matrix is not so good for transitions with 
much smaller matrix elements such as 4s — ?> Qpj. 

The ratio of the line strengths for the 4s — >■ Apj tran- 
sition would be exactly 2.0 in a non-relativistic model. 
Experiment and theory both indicate that the difference 
of the 4s — Ap transition ratio from 2.0 is very small. 
The current calculation and the MBPT-SD calculation 
indicate that the matrix element ratio is slight smaller 
than 2.0. However, the ratio differs substantially from 
2.0 for the 4s — >■ 5pj and 4s 6pj transitions. The 
main cause for the deviation of the ratio from 2.0 is the 
slightly different wavefunction for the spin-orbit doublet 
arising from the slightly different energies [23|. We have 
done non-relativistic calculations and have been able to 
reproduce the ratios given in Table IIIII by simply tuning 
the polarization potential to separately give the experi- 
mental binding energies of each spin-orbit doublet. 



C. Polarizabilities and tune-out wavelengths 

The computations of the static polarizabilities utilized 
the RCICP matrix elements but with the excitation en- 
ergies for the 4pj, 5pj and 6pj set to be those of exper- 
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TABLE III: Comparison of reduced dipole matrix elements 
(a.u.) for the principal transitions of potassium with experi- 
mental values and other calculations. 

Transition RCICP MBPT-SD CICP Experiment 



[13,i8,i9, 23] 


4s - 4pi/2 


4.1030 


4.098 


4.1017 


4.102(5) [24] 








4.106(4) [25] 


4s - 4p3/2 


5.8016 


5.794 


5.8006 


5.800(8) [24] 








5.807(7) [25] 


4s - 5pi/2 


0.2634 


0.271(5) 


0.2696 




4s - 5p3/2 


0.3886 


0.398(8) 


0.3812 




4s - 6pi 10 


0.0756 


0.084 


0.0791 




4s - 6p3/2 


0.1162 


0.128 


0.1118 




4pi/2 -5s 


3.8879 


3.855(1) 


3.9058 




4p3/2 -5s 


5.5384 


5.54(1) 


5.5236 




5s - 5pi/2 


9.4967 


9.49(3) 


9.4918 




5s - 5p3/2 


13.410 


13.40(4) 


13.423 




5pi/2 -6s 


8.7766 


8.79(2) 


8.8088 




5p3/2 -6s 


12.490 


12.50(2) 


12.458 




4pi/2 -3^3/2 


7.9662 


7.97(3) 


7.9812 


7.979(35) [13] 


4p3/2 -3^3/2 


3.5731 


3.57(1) 


3.5693 


3.578(16) [13] 


4p3/2 -3^5/2 


10.719 


10.72(3) 


10.708 10.734(47) fl3] 


4pi/2 -4^3/2 


0.1403 


0.1121(8) 


0.1246 




4p3/2 -4^3/2 


0.0529 


0.0400(1) 


0.0557 




4p3/2 -4^5/2 


0.1553 


0.1170(4) 


0.1672 




3^3/2 -5pi/2 


7.1687 


7.2(1) 


7.1476 




3^3/2 -5p3/2 


3.1934 


3.19(9) 


3.1965 




3^5/2 -5p3/2 


9.5743 


9.6(1) 


9.5895 




5pi/2 -4^3/2 


17.040 


17.04(6) 


17.076 




5p3/2 -4^3/2 


7.6432 


7.64(3) 


7.6367 




5p3/2 -4^5/2 


22.932 


22.93(8) 


22.910 
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1.99935 1.9987 [5] 2.0 



2.17787 



2.35799 



2.17964 



2.31894 



2.0 



2.0 



2.0005(40) [5] 
1.999(4) [24] 
2.002(3) [25] 
2.01 [26] 

2.15 [26] 
2.28 [261 



iment. The dynamic polarizability is defined 

/on 



(6) 



where /on is the oscillator strength for the dipole transi- 
tion. For low frequencies, the dynamic polarizability can 
be expanded 



i{uj) = a(0) + 5(-4)cj2 + 



(7) 



where a(0) is the static dipole polarizability calculated 
at = and 5 (—4) is calculated using the oscillator 
strength sum rule, 



5(-4) = E% 



(8) 



Polarizabilities for the potassium ground state from dif- 
ferent sources are listed in Table ITVl The present RCICP 



TABLE IV: Static dipole polarizabilities (in a.u.) for potas- 
sium. A short description of the details behind some of the 
polarizabilities can be found in Ref. 





Ql 


Q2 (10^) 


Theory: Present RCICP 


290.1 


5.000 


Theory: CICP [Tj 


290.0 


5.005 


Theory: MBPT-SD [19] 


289.3 


5.018 


Theory: RCCSD [30J 


301.28 


5.018 


Expt. Ey.H [31] 


293(6) 




Expt: Interferometer ^8] 


290.8(1.4) 




Hybrid: Sum rule [29] 


290.2(8) 





TABLE V: Pseudo-spectral oscillator strength distribution for 
the potassium core. Energies are given in a.u.. 



133.6890020 
14.6459330 
11.6752580 
1.90477720 
1.11041710 



2.0 
2.0 
6.0 
2.0 

6.0 



calculation gave essentially the same polarizability, 290.1 
flg as its non-relativistic CICP predecessor. The non- 
relativistic CICP model had its valence energies tuned to 
experimental energies (in this case, the (2J-I-1) weighted 
energy of any spin-orbit doublet), just like the present 
calculation. The As—Apj radial matrix elements are dom- 
inated by long range form of the wave function, and that 
is practically the same for the CICP and the relativis- 
tic RCICP calculations. The present RCICP calculation 
agrees with experiment (28l . [29l ] within the experimental 
uncertainties. 

The relativistic coupled cluster calculation (RCCSD) 
[30| gives a dipole polarizability that is about 3% larger 
than experiment. 

The core polarizabilit y is given by a pseudo-oscillator 
strength distribution 0,112, [1^. The distribution is tab- 
ulated in Table |Vl The distribution is derived from the 
single particle energies of the Hartree-Fock core. Each 
separate (n, £) level is identified with one transition with 
a pseudo-oscillator strength equal to the number of elec- 
trons in the shell. The excitation energy is set by adding 
a constant to the Koopman energies and tuning the 
constant until the core polarizability from the oscillator 
strength sum rule is equal to the know core polarizability. 
In the present case, the K+ core polarizability was set to 
5.47 a.u. 34, 35]. 

The tune-out wavelengths all tend to be close to the 
wavelengths for excitation of the npj excited states. 
There are two scenarios that lead to tune out wave- 
lengths. In the first, the tune-out wavelength occurs 
in the middle of an npj spin-orbit doublet. The wave- 
length will be shorter than the transition wavelength to 
the npi/2 state and longer than the wavelength to the 
np3/2 state. When this occurs, the dynamic polarizabili- 
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ties of the npi/2 and states will have the opposite 

sign and this will lead to a zero in the total dynamic 
polarizability occurring for all spin-orbit doublets. The 
second scenario leading to a tune-out wavelength occurs 
when the wavelength is shorter than that for excitation 
of the Apj states. When this occurs, the contribution to 
the dynamic polarizability from the Apj states become 
negative. This leads to a series of tune-out wavelengths 
occurring just below the excitation energies of the 5pi/2, 
6pi/2, 7pi/2, . . .states. 

Simplified expressions can be used to describe the dy- 
namic polarizabilities on the vicinity of the tune-out 
wavelengths. The first tune-out wavelength occurs when 
the wavelength lies between the 4px/2 ^nd 4^3/2 reso- 
nant wavelengths. The dynamic polarizability here can 
be written, 



4pi/2 



4p3/2 



(Ai?2 



4pi/2 



4p3/2 



■ + arcm(w) , 

(9) 

where a^cmi'-^) is remainder part of dynamic dipole polar- 
izability. The energy difference, AE4p^.^ can be param- 



eterized as AK 



4p3/2 



AE4p^^^{l + S). Parameterizing 
the line strength, S (the square of the reduced matrix 
element), as 

S{As ^ 4p3/2) = S{As ^ 4pi/2)(2 + R) (10) 



leads to 



ai(w) 



/4pi 



/2 



/4p,,,(2 + i?)(l + (5) 



(AEj^^^^-oj^) [{AE4p,,,ni + dr-u^] 

+ arom('^) • (11) 

The dipole oscillator strength f^p-^^^ obtained by multi- 
plying the reduced matrix element, with the experimen- 
tal 4s Ap>Yj2 energy difference. The value of R is sim- 
ply the ratio of computed line strength coming from the 
RCICP calculations. It can be seen from Table IVTl that 
the remainder polarizability, Q!rcm(w) only makes a small 
contribution to the total polarizability. The remainder 
polarizability varies relatively slowly with wavelength in 
the vicinity of the tune-out wavelength. 

Table I VIII illustrates the variation in the tune-out 
wavelength with respect to variations in R and arcm(<^)- 
The contributions to the polarizability from the 4p2/2 
and 4p3/2 transitions are 5,000 times larger than those 
from every other transition. A change in arem(w) of 1.0 
a.u leads to the tune-out wavelength changing by 0.00002 
nm. The tune-out wavelength is much more sensitive to 
variations in R. A change in R to —0.005 leads to the 
tune-out wavelength changing by 0.0016 nm. A value of 
R = -0.005 is 3 times larger than the MBPT-SD value 
and is seven times larger than the RCICP value of R. 

A different parameterization should be used in the 



vicinity of the of the excited states with n > 4. 

/npi/2 



ai(w) = a4p('^) + arcm('^) + 



(A R2 - u;2 



/„p,,,(2 + j?)(l + ^) 

[(A£;„p^,J2(l + 5)2_^2] 



(12) 



Here the polarizability arising from the 4s — )> 4pj transi- 
tions is retained as a separate term since it is much larger 
than the remainder. 

The tune-out wavelengths in the vicinity of the 5pj 
levels illustrate clearly how the arem(w) and R are of 
different importance depending on whether the tune-out 
energy is located between the 5pi/2 and ^p^ii levels or 
before the 5pi/2 level. 

Table IVII shows that the tunc-out wavelength at en- 
ergies below the 5pi/2 excitation threshold is caused by 
the cancellation of the hpj and core contributions with 
those of coming from the 4pj levels. This tune-out wave- 
length is effectively determined by ratio of the 4p and 
hp oscillator strengths. The contribution of the core is 
small in absolute terms, and a 5% uncertainty in the 
core polarizability will have a small effect on the tune- 
out wavelength. The tune-out wavelength is predomi- 
nantly determined by the relative size of the 4s Apj 
and 4s — > 5pj matrix elements. This can be seen from 
Table IVIII Increasing the 4s 5pj matrix elements by 
3% (roughly the difference with the MBPT-SD matrix 
elements) leads to the tune-out wavelength increasing by 
0.077 nm. Changing the value of R, by 0.02 leads to the 
tune-out wavelength changing by 0.007 nm. It should 
be noted that our definition of R does imply an overall 
increase in the total oscillator strength to the 5pj states. 

The tune-out wavelength at the energies between the 
5pi/2 and 5p3/2 levels does have some dependence on 
arem siuce it now incorporates the contribution from the 
4s — ^ 4pj transitions. A 5% change in the matrix element 
leads to a change in the tune-out wavelength of 0.0024 
nm, this is 60 times smaller than the effect on the tune- 
out wavelength at 405.9173 nm. The sensitivity to a 
change in the value of R by 0.02 was only 0.0007 nm. In 
relative terms, the tune-out wavelength is more sensitive 
to the value of R than |(5pj||il'||4s)| in the 5pi/2 to 5p3/2 
energy gap, than it is in the energy region before 5pi/2 
excitation. 

Table IVIII also gives the tune-out wavelengths in the 
vicinity of the 6pj levels. Once again the tune-out wave- 
lengths are very sensitive to the absolute size of the 
|(6pj||-D||4s)| transition matrix element. A 10% change 
in the matrix element (the difference between the RCICP 
and MBPT-SD calculations) leads to a change of 0.0234 
nm in the tune-out wavelength just below the 6pi/2 
threshold. The considerations that determine the val- 
ues of the tune-out wavelengths in the vicinity of the 
bp J states also apply to the tune-out wavelengths in the 
vicinity of the 6pj states. 

Finally, it is noted that the tune-out wavelengths were 
also evaluated using the MBPT-SD matrix elements in 
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TABLE VI: Breakdown of contributions to the potassium ground state polarizability at different wavelengths. 



A (nm) 


oo 


768.97075 


405.9173 


404.7217 


344.9099 


344.7861 


cj (a.u.) 





0.059252386 


0.11224787 


0.11257945 


0.13210218 


0.13214964 


4s - 4pi/2 


94.8454 


-32032.4796 


-36.4876 


-36.1911 


-23.7987 


-23.7774 


4s - 4p3/2 


188.7902 


32025.7787 


-73.5278 


-72.9282 


-47.9006 


-47.8575 


4s - 5pi/2 


0.2054 


0.2842 


38.5906 


—364.9459 


—0.5438 


—0.5424 


4s - 5p3/2 


0.4469 


0.6180 


65.4323 


468.0690 


-1.1901 


-1.1870 


4s - 6pi/2 


0.0144 


0.0181 


0.0519 


0.0527 


33.7885 


-49.5664 


4s - 6p3/2 


0.0340 


0.0426 


0.1221 


0.1240 


33.7302 


117.0156 


Remainder Valence 


0.2426 


0.2528 


0.2957 


0.2963 


0.3712 


0.3717 




5.4708 


5.4852 


5.5229 


5.5232 


5.5433 


5.5434 


Total 


290.050 


















TABLE VII: Values of the tune-out wavelength for the K 
atom. The experimental transition wavelengths are taken 
from [3^ . [37I I . The first entry lists the tune-out wavelength 
as computed with RCICP matrix elements. The other en- 
tries exhibit the changes to the tune-out wavelengths when 
changes are made to matrix elements underlying the oscilla- 
tor strength sums. 



Resonance 



oj (a.u.) 



A (nm) 



A£4s-^4pi/2 0.059164859 770.10836 

0.059427807 766.70089 
R = -0.0006476, Qrom(t^) = 6.701 0.059252386 768.97077 
arem{i^) = 5.701 0.059252387 768.97075 

R = -0.005 0.059252513 768.96912 

MBPT-SD [4] 0.0592524(2) 768.971(3) 

Experiment [5] 0.0592523(1) 768.9712(15) 



Ai54s-5p3/2 

R = 0.17787, arcm(cj) 
R = 0.15787 
K5pj||D||4s)| X 1.03 
MBPT-SD [4] 



-104.023 



11254778 

11263324 

11224787 

1122499 

1122120 

11223(1) 



404.8356 
404.5285 
405.9173 

405.9100 
405.9943 
405.98(4) 



R = 0.17787, arcn,{ui) = -103.300 0.11257946 404.7217 

7£ = 0.15787, 0.11257967 404.7210 

K5pj||L>||4s)| X 1.03 0.11257898 404.7228 

MBPT-SD [4] 0.11258(1) 404.72(4) 

AB4S-6P1/2 0.13213040 344.8363 

AE4s-6py2 0.13216885 344.7359 

R = 0.35799, arcmioj) = -67.5188 0.13210218 344.9099 

R = 0.30799, 0.13210256 344.9089 

|{6p,/||D||4s)| X 1.10 0.13209360 344.9323 

MBPT-SD [4] 0.1320933(4) 344.933(1) 



R = 0.35799, Qfrcm 

R = 0.30799, 
(6pj||D||4s)| X 1.10 



-67.4492 0. 

0. 
0. 



13214964 
13214992 
13214825 



344.7861 

344.7853 
344.7897 



Table IIIII but with other aspects of the calculation taken 
from the RCICP calculation. The resulting tune-out 
wavelengths were identical to the MBPT-SD tune-out 
wavelengths in Table IVlIl to all quoted digits. 



TABLE VIII: Tabulation of atomic parameters, with esti- 
mated uncertainties that can be used reduced dipole matrix 
elements (a.u.) for the principal transitions of potassium with 
experimental values and other calculations. 



Parameter 



RCICP 



4s - 4pi/2 
4s - 4p3/2 
4s - 5pi/2 
4s - 5p3/2 
4s - 6pi/2 
4s - 6p3/2 

dcorc — valence 

Remainder Valence 
emainder Valence 

(-4) 



4.1030(7) 
5.8016(10) 
0.2634(8) 
0.3886(100) 
0.0756(100) 
0.1162(14) 
5.4708(1000) 
4.10(80) 
-0.13 
0.243(111) 
2.67(134) 



D. Uncertainties 

Part of this manuscript is focussed on the prediction 
of the tune-out wavelengths, but another and possibly 
more important part concerns the extraction of useful 
atomic structure information from an experimental value 
of the tune-out wavelength. Knowledge of the npj tune- 
out wavelengths permits the determination of the npj 
oscillator strengths to a high degree of precision. 

The most important atomic parameters that con- 
tribute to the long wavelength dynamic polarizability 
are listed in Table IVIIII These parameters are derived 
from the RCICP calculations and uncertainties are esti- 
mated by examination of the difference with experiment 
or MBPT-SD calculations. 

The determination of the 4pi/2 : ^P3/2 line strength 
ratio is only weakly dependent on the value of the non 
4p terms in the dynamic polarizability since these terms 
are small. 

Knowledge of the npj tune-out wavelengths permits 
the determination of the npj oscillator strengths for 
n > 5 to a high degree of precision. The polarizability 
becomes zero when the contributions to the polarizabil- 
ity from the npj levels and the remainder cancel exactly. 
The biggest terms in the remainder are a^pj and acoro 
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which together constitute 99% of the remainder. More- 
over, both of these terms are known with a reasonable 
degree of precision. 

The uncertainties in the experimental 4s Apj line 
strengths do not exceed 0.3% and the RCICP transition 
matrix element lies between two experimental estimates 
[M HI- There is 2% variation between the RCICP and 
MBPT-SD estimates of the 4s — )■ 5pj matrix elements 
and a 10% variation between 4s 6pj matrix elements. 
Taken together, 4s — >■ 5pj and 4s — ?• 6pj would con- 
tribute less than 0.01 a.u. to the uncertainty of the total 
polarizability at uj — a.u. 

The uncertainty in the core polarizability of 5.47 a.u. 
itself is stated to be about 2% [H]. This polarizability 
was based on the binding energies of the 4/, 6/ and 9/ 
levels of potassium using spectral data from 1955 j3^]. 
There is scope for an improvement in the precision of the 
core polarizability and this could be easily accomplished 
by spectroscopic experiments that measured the energies 
of the ng levels. The construction of pseudo-oscillator 
strength distribution for the core permits the energy 
variation of acore to be incorporated into the calculation. 
The uncertainty of 20% was estimated by using the same 
procedure to construct the pseudo-oscillator strength dis- 
tribution to argon and making reference to a hig hly ac- 
curate pseudo-oscillator strength distribution ^3^. 

The core- valence term, acoro-vaicncc is a term that com- 
pensates for Pauli-Principle violating excitations from 
the core to the valence 4s orbital. This value is sourced 
from an MBPT-SD calculation [23] since it is not incor- 
porated in the RCICP calculation. No uncertainty has 
been assigned to this contribution to the polarizability. 

The valence remainder term contains contributions 
from highly excited discrete transitions as well as contri- 
butions from the continuum. The present value is 0.243 
a.u. a.t u ~ 0. This is more than twice the size of a 
MBPT-SD estimate of 0.07 a.u. However, much of 

the MBPT-SD valence remainder is computed in the DF 
approximation. We have performed calculations in the 
DF approximation, and the DF oscillator strengths em- 
bedded in the continuum beyond the cooper minimum at 
0.010 a.u. dl] are typically a factor of 3 and 4 smaller 
than the RCICP oscillator strengths at those energies 
Nevertheless, the uncertainty in this term has been con- 
servatively assessed at 50%. Comparisons of the RCICP 
and MBPT-SD oscillator strength distributions in the 
continuum, and further comparisons with experimental 
photo-ionization cross sections would be helpful in re- 
fining the estimate and uncertainty of this rather small 
term. 

The parameter set and error budget in Table I Villi can 
be utilized to help convert tune-out wavelengths into os- 
cillator strengths. There is room for improvement in 
the parameter set. Measurements of the tune-out wave- 
lengths near the 5pj excitation will result in better esti- 
mates of the 5pj matrix elements. Parameters obtained 



from theory do not have to be exclusively obtained from 
a single calculation, for example some information from 
Table|VllI]might be best obtained from a MBPT-SD cal- 
culation while others, e.g. the valence remainder might 
be best estimated from the RCICP calculation or some 
RCICP/MBPT-SD hybrid. 



IV. CONCLUSION 

The five lowest tune-out wavelengths for the potassium 
atom are computed by a relativistic structure model. A 
detailed analysis is performed regarding the contribution 
that the different terms make to the polarizability. The 
results illustrate the dependence of the tune-out wave- 
lengths on a relatively small number of atomic parame- 
ters. 

The lowest energy tune-out wavelength is primarily 
determined by the ratio of the line strengths for the 
4s — Apj transitions. The present calculation, and the 
MBPT-SD calculation are in agreement with existing ex- 
perimental data [53. The precision of the experiment 
would need to improve by an order of magnitude to pro- 
vide a stringent test of the Apj states line strength ratio. 
However, Holmgren et. al. |5| suggest that it might be 
possible to improve the precision by up to 3 orders of 
magnitude. 

The tune-out wavelengths near the 5pj excited states 
are most sensitive to the ratio of the 4s Apj and 
4s 5pj matrix elements. The remainder term in- 
corporating all transitions except for the 4s bpj is 
dominated by the 4s Apj and core polarizabilities 
which comprise 99.5% of the remainder polarizability. 
The tune-out wavelengths here provide a means to deter- 
mine the 5s 5pj oscillator strengths to high precision. 
Measurement of the tune-out wavelength to a precision 
of 0.01 nm would lead to oscillator strengths with a pre- 
cision better than 1%. 

Holmgren et. al. 0] suggested that measurements of 
the tune-out wavelengths near the 5pi/2 excited states 
could be used to determine the core polarizability. We do 
not agree with this statement. The remainder terms near 
the 5pj excitation are dominated by a4p{uj) and acoic('^)- 
The uncertainty in a4p{uj) at 405 nm would be about 0.3 
a.u. and this uncertainty limits the precision with which 
the core polarizability could be measured. The preferred 
approach is to treat the a4p{uj) and Q;coro(w) polarizabil- 
ities as known quantities with relatively small uncertain- 
ties and use measurements of the tune-out wavelength to 
extract precision values of the 5pj oscillator strengths. 

This research was supported by the Australian Re- 
search Council Discovery Project DP-1092620. We thank 
Dr M S Safronova for helpful communications regarding 
her matrix elements. 
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